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Abstract. A ROSAT HRI observation of the cluster MS- 
0839.9+2938 at z=0.194 is presented. It confirms the earlier 
suggestion, based on the detection of extended Ha emission, 
that the inner regions of this cluster are dominated by a cool- 
ing flow. The surface brightness distribution within the cooling 
radius shows structures which can be interpreted as evidence 
of a highly inhomogeneous space distribution of the cooling 
process. We note that the brightness at the barycentre of its 
distribution, which falls right on top of the central giant ellip- 
tical galaxy, is lower than the peaks in the structures around: 
we suggest that this situation is likely to arise as a consequence 
of photoelectric absorption by cold gas within the cooling flow, 
with an equivalent column density in the order of 5x 10 21 cm -2 
within ~10" from the centre. 

Key words: Galaxies: clustering; cooling flows; Galaxies: clus- 
ters: individual: MS0839+2938 
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The cluster MS0839+2938 was discovered as an extended X-ray 
source in the EINSTEIN Extended Medium Sensitivity Survey 
(Gioia et al. 1987, 1990). Subsequent optical and radio obser- 
vations (Nesci et al. 1989) showed that the cluster is domina- 
ted by a giant elliptical, probably a cD type galaxy, which is 
characterized by an intense and extended Ha emission, and a 
compact, relatively weak radio emission. The redshift of this 
galaxy is z=0.193, and the cluster galaxy velocity dispersion 
along the line of sight, estimated on the basis of six redshifts, 
is 1500±700 km s" 1 . At the distance of 1200 Mpc (H =50 km 
s -1 Mpc -1 , qo = 0.5), the X-ray luminosity in the 0.3-3.5 keV 
band was estimated to be about 4xl0 44 erg s -1 . 

Since a spatially extended optical line emission is present in 
several nearby cooling flow clusters, MS0839+2938 appeared 
to be one of the most distant candidate cooling flow cluster 
known at the epoch of the ROSAT launch. Strictly speaking, 
the identification of a cooling flow would require the detection 
in the X-ray emitting gas of a much lower temperature in the 
innermost regions than in the whole cluster. However, when 
this goal cannot be directly achieved, a cooling flow regime is 
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generally regarded as proven if the central density, as derived 
from the cluster surface brightness distribution adopting the 
dominant cluster temperature, implies a cooling time definitely 
shorter than the Hubble time (see e.g. Arnaud 1988). The aim 
of our ROSAT observation was therefore to derive the density 
profile of MS0839+2938 and verify the existence of the cooling 
flow. Since the scale at the cluster distance is 4.13 kpc arcsec -1 , 
in order to resolve the emission within the core radius we used 
the HRI instrument. This choice, however, implied the loss of 
temperature information. A detailed description of the satellite 
and its instruments can be found in the "ROSAT A02 Call for 
Proposals" (MPE, 1991) 



2. Data analysis 

The pointed observation of the cluster MS0839 +2938 (image 
number WG800159H.N2) lasted 18879 live seconds and was 
performed on April 20-22, 1992. Besides the extended image 
of the cluster at the center of the field of view, a point source 
was also detected 4.7' to the SW (32 net cts), corresponding 
to the star HD74076; this allowed us to check the zero point 
of the sky coordinates of our image and to overlay the X-ray 
onto the optical image of the cluster with the accuracy of ~2" . 
We determined the background level by measuring the 
counts in 32 adjacent squares with sides of 260" all around 
the cluster in the central, 1600" side square of the frame. The 
average count rate is 1.150xl0 -6 cts arcsec -2 s -1 (0.0217 
cts arcsec -2 ), very near to the pre-flight expected value of 
1.1 x 10 -6 and to the value (1.25 x 10 -6 ) found by Sarazin et al. 
(1992a) in their ROSAT HRI observation of 2A 0335+098. The 
standard deviation of the 32 measurements is fully consistent 
with the Poisson statistics. We also measured the background 
level in a ring 200"-280" around the cluster center finding the 
same result. 



3. X-ray morphology of the cluster 

In order to evaluate the maximum extent out to which the 
cluster emission is detected, after background subtraction an 
isophotal map was obtained by smoothing the original data, 
binned in 1" x 1" pixels, through a Gaussian filter with a = 



32". The map is shown in Fig. 1, where the lowest contour 
corresponds to four times the background fluctuations: at this 
brightness level (0.005 cts arcsec~ 2 ) the emission is elongated 
in the N-S direction, as the dominant galaxy, with a major axis 
of ~5 arcmin (~1200 kpc). An approximately circular type 
of symmetry is evident in the less smoothed Fig. 2, obtained 
through a Gaussian filter with a = 8" , which shows the emis- 
sion to be centrally peaked, with the peak coinciding with the 
dominant galaxy, thus leading to a morphological classification 
as an XD cluster, according to the scheme proposed by Jones 
and Forman (1984). 



Fig. 1. The X-ray map at low resolution after background subtrac- 
tion. The isophotes arc at levels of 0.0051, 0.0076, 0.0102, 0.0127, 
0.0191, 0.0254, 0.0317, 0.0381, 0.0508 net cts arcsec~ 2 . 



The isophotal map at high resolution (unbinned and smoo- 
thed with a Gaussian filter with a =2") of the very central 
region is shown in Fig. 3 overlaid onto an optical (R) image of 
the cluster (Nesci et al. 1989). This map shows that the bright- 
ness distribution in the central region is not smooth, but that it 
contains structures with sizes in the range of a few tens of kpc. 
The brightness contrast and the length scale of these struc- 
tures are consistent, taking into account the blurring due the 
difference in distance, with those discovered with the ROSAT 
HRI in two relatively nearby cooling flow clusters, 2A0335+096 
(z=0.035) and A2029 (z=0.0767) by Sarazin et al. (1992a, b). 
Given the limited count statistics, the details in the observed 
structures, such as their exact brightness excess, can be deter- 
mined only with marginal significance, except for the bright 
and resolved blob located about 12" to the South of the giant 
elliptical: the peak excess brightness of this blob is approxi- 
mately equal to that of the surrounding diffuse emission, and 
the size is about 6", or 24 kpc, across. 

Another important morphological aspect is that the bright- 
ness, at the position of the giant elliptical galaxy, is lower than 
the peaks in the surrounding structures, while the radio posi- 
tion of the galaxy centre (8 h 42 m 55.9 s ; +29°27'27" (Eq. 2000)) 
coincides within about 2" with the barycentre of the overall 



Fig. 2. The isophotal map at intermediate resolution: the isophotes 
arc at levels of 0.024, 0.044, 0.064, 0.084, 0.104, 0.144, 0.184 and 
0.264 net cts arcsec -2 . 



Fig. 3. Map at the highest resolution of the inner cluster region 
overlaid onto to the optical image of the cluster. The isophotal con- 
tours arc 0.1383, 0.2183, 0.2983, 0.3783, 0.4583, 0.5383 and 0.6183 
net cts arcsec — 2 . 



X-ray emission. We shall return on these aspects in Sections 4 
and 5. 

As a side issue, we note that no X-rays were detected from 
the emission line galaxy (G27 in Nesci et al. 1989) with a very 
blue colour (V-R=0.42), located 68" West of the dominant 
cluster galaxy and at a similar redshift (z=0.186). 



4. Physical parameters 

The radial profile of the brightness distribution, azimuthally 
averaged in steps of 10", is given in Fig. 4. The profile be- 
yond 50" (^200 kpc) can be well fitted with a power law of 
slope s=— 2.12±0.42 at the 95% confidence level. This slope 
corresponds to a value for the (3 parameter of the hydro- 
static isothermal model (Cavaliere and Fusco-Femiano 1976) of 
0.52±0.07, similar to those found by Sarazin et al. (1992a, b) 
using the ROSAT HRI for 2A 0335+096 (0.55) and A2029 
(0.61). 



Fig. 4. The profile of the surface brightness 5" azimuthally averaged 
in rings of 10" , as a function of radius r. The best fitting power-law 
profile beyond of 50" is shown as a solid line. Error bars represent 
the 68% confidence level in Poisson statistics. 



The conversion from observed counts to physical flux re- 
quires an assumption for the X-ray emission spectrum of the 
cluster and an estimate of the absorption along the line of 
sight. For the plasma thermal spectrum we adopted the code 
by Mewe et al. (1985) with 0.3 solar abundance, a typical value 
for intracluster gas (Edge et al. 1991a), and folded it with the 
effective area of the ROSAT HRI. A correction was also applied 
to take into account the energy band shift due to the cluster 
redshift (k-correction) , but this correction is very small with 
respect to that due to the galactic absorption. For the latter a 
hydrogen column density of log(Nn)=20.6 was adopted, as in 
Nesciet al. (1989). 

The gas temperature T can be evaluated using the empirical 
correlation (Edge and Stewart 1991b) between T and the ra- 
dial velocity dispersion of the cluster galaxies. However, given 
the small number (6) of redshifts available, the dispersion esti- 
mate is subject to a large uncertainty related to projection, and 
possibly even contamination effects. We therefore preferred to 
adopt two rather different values of kT (3 and 7 keV) and to 
check 'a posteriori' whether the derived bolometric cluster lu- 
minosity fits in the empirical luminosity-temperature relation 
(Edge and Stewart 1991a). 



The absorption corrected flux within 120" from the cluster 
centre, in the band 0.1-2.4 keV, is (2.45±0.11)xl0~ 12 erg cm~ 2 
s" 1 for kT=7 keV, and (2.50±0.11) x 10~ 12 for kT=3 keV: these 
values are consistent with the estimate based on the EINSTEIN 
Observatory observation within the same radius, (2.23±0.24) 
x 10~ 12 erg cm~ 2 s" 1 rescaled to the ROSAT band. Within 
a radius of 200", which encircles the lowest contour of Fig. 
1, and contains 948 net counts, the total flux is (2.83±0.14) 
xl0~ 12 ergem^s- 1 for kT=7 keV, and 2.86xl0~ 12 for kT=3 
keV. This corresponds to a luminosity in the ROSAT band of 
4.8xl0 44 erg s _1 for the temperature range considered; the 
bolometric luminosity then is 1.25xl0 45 erg s _1 if kT=7 keV, 
or 7.5 xlO 44 erg s" 1 if kT=3 keV. 

Edge and Stewart (1991a, their Fig. 11), showed that there 
exists a correlation between the temperature and the bolome- 
tric luminosity, whose normalization however depends on the 
value of the central density. As we shall presently see, the cen- 
tral density derived from the X-ray data is rather insensitive to 
which one of the two temperature values is adopted, and our 
cluster would belong to their subsample of objects with cen- 
tral density > 9xl0~ 3 cm~ 3 . With respect to their best fitting 
correlation line in the kT-Lt, ! plane, this cluster falls to the 
right (its luminosity is too low) if kT=7 keV, to the left (its 
luminosity is too high) if kT=3 keV. The actual cluster tem- 
perature is therefore likely to be in this range. In the following 
we shall derive the relevant parameters for both cases. 

The cluster one dimensional velocity dispersion expected 
from the derived value of (3 is ~800 km s _1 for kT=7 keV and 
~520 km s _1 for kT=3 keV. The estimate based on six redshifts 
and reported in Section 1 is marginally consistent with the first 
value, but it could be overestimated due to the contamination 
by a non-member object. After removing the most discrepant 
redshift we obtain a dispersion of 850 (±400) km s -1 , quite 
consistent with either of the two values derived above. 

From the observed surface brightness profile, assuming 
spherical symmetry and a constant temperature, it is possible 
to derive the density profile with a deprojection technique (see 
e.g. Arnaud 1988). To overcome numerical problems arising 
from the poor statistics on the points at large radii, we de- 
projected an analytical fit of the brightness profile. The most 
relevant resulting parameters, namely the central density no, 
the central cooling time t co , the cooling radius r c (defined as 
the radius where the cooling time equals the Hubble time mi- 
nus the look back time, that is 10 10 years), are given in Table 
1 for the two temperatures; the density profile is illustrated in 
Fig. 5. 



Table 1: Physical parameters 




kT no t co r c 
(keV) (cm -3 ) (Gyrs) (kpc) 


M 9as (<0.5 Mpc) 
(M ) 


7 0.0170 3.96 70 
3 0.0156 2.57 110 


l.lxlO 13 
1.2xl0 13 



From this table it can be seen that the derived central den- 
sity depends very weakly on the adopted temperature. Further- 
more, no matter which one of the two temperature values is 
adopted, the central cooling time of MS0839+2938 is definitely 
lower than the Hubble time, so this cluster certainly contains 
a cooling flow. The mass flow rate within the cooling radius 
can be estimated from the relation M — 0.4/im_srLj, o ;/(fcT), 
where jj, is the mean molecular gas weight, Lboi is the bolome- 
tric X-ray luminosity within the cooling radius and nin is the 



for kT=7 keV and 5. Discussion 



hydrogen mass. We derive M—130 Mq y 
M=280 M y" 1 for kT=3 keV. 

It is noteworthy that the map in Fig. 3 is dominated by 
emission from the cooling flow region, hence the structures vi- 
sible in the map lie within the cooling radius and are likely 
to be the result of inhomogeneities where the cooling process 
goes on more rapidly. By assuming that these structures are 
in pressure equilibrium with the surrounding medium, one can 
estimate their temperature and therefore their proper cooling 
time. Let us for instance consider the bright blob to the South 
of the central galaxy, and assume that it has a spherical con- 
figuration with radius 12 kpc and a filling factor equal one. 
By adopting a distance from the centre equal to the projected 
distance, from Fig. 5 one obtains the pressure term nT; for an 
external temperature of 7 keV, it turns out that the observed 
excess counts in the blob require that the internal temperature 
be about 3.5 keV, and the internal density therefore be twice 
the external value. The cooling time in the blob is then about 
2x 10 9 years, a factor of four shorter than that of the surround- 
ing gas. The temperature would be lower, and the cooling time 
shorter, if the filling factor were less than unity. 



Fig. 5. The proton density profile n p for kT=7 keV. The position of 
the cooling radius r c is marked. For kT=3 keV the density is lower 
by 0.04 dcx. 



The ROSAT observation described in this paper demonstrates 
that the inner regions of the distant cluster MS0839. 9+2938 are 
dominated by a cooling flow, and that, very much like what 
is observed in two relatively nearby clusters (Sarazin et al. 
1992a, b), the space distribution of the cooling process is highly 
inhomogencous. 

Here we would like to address the question whether photo- 
electric absorption by cool gas within the cooling flow might be 
affecting the observed surface brightness distribution. White et 
al. (1992) present spectral evidence, based on EINSTEIN Ob- 
servatory SSS data, that cold, X-ray absorbing gas is often 
present in cooling flows, with an equivalent N_s up to a few 
times 10 21 cm~ 2 : such a column density would significantly af- 
fect the ROSAT HRI counts. They argue on physical grounds 
that this gas is likely to be in the form of dense and small 
clouds. If these clouds are regarded as remnants of a cooling 
process which went on for 10 10 years, even if the space distri- 
bution of their place of origin were highly inhomogeneous, it is 
quite conceivable that their present, cumulative distribution is 
a regular function of the radial distance from the centre. There- 
fore the associated absorption can hardly be responsible for the 
structures in the surface brightness observed within the cooling 
radius. On the other hand, as it was noted in Section 3, the 
brightness level at the barycentre of its distribution, where the 
central giant elliptical is located, is lower than the peaks in the 
surrounding structures. We suggest that this situation could 
be due to photoelectric absorption. In fact, a column density 
of 5xl0 21 cm~ 2 would reduce the brightness in the ROSAT 
band by a factor of about 2 for kT >2 keV. A simple calcula- 
tion shows that, if M had remained constant over 10 10 years, 
and all the cooled gas went into X-ray absorbing clouds, the 
Nh that one obtains if the clouds were uniformly distributed 
within the cooling radius, and with a covering factor of unity, is 
7xl0 21 cm~ 2 . This value is a lower limit to that which would 
be obtained if the clouds were more realistically assumed to 
concentrate towards the centre. We conclude that a spectral 
imaging of the cluster will most probably reveal the presence 
of a strong absorption within 10" or so from the centre, a goal 
achievable in the near future with the angular resolution of the 
AXAF satellite. 
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